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Abstract 

We investigate the dependence of the nucleon-nucleon force in the deuteron 
system on the values of coupling strengths at high energy, which will in 
general depend on the geometry of extra dimensions. The stability of deu¬ 
terium at all times after nucleosynthesis sets a bound on the time variation 
of the ratio of the QCD confinement scale to light quark masses. We hnd 
the relation between this ratio, which is exponentially sensitive to high- 
energy couplings, and fundamental parameters, in various classes of unified 
theory. Model-dependent effects in the Higgs and fermion mass sector may 
dominate even over the strong dependence of the QCD scale A. The bind¬ 
ing energy of the deuteron also has an important effect on nucleosynthesis: 
we estimate the resulting bounds on variation of couplings. 


1 Introduction 

In many models of particle physics, the universe is assumed to have more than 
4 dimensions. The extra dimensions are either compactihed to such a small size 
that we cannot (currently) probe them experimentally PP or possess metrics with 
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a nontrivial dependence on the transverse directions such that we can only de¬ 
tect the gravitational influence of our familiar 4 dimensions jS!- Cosmological 
solutions of the held equations of these theories often involve time evolution of 
the higher dimensions; the value of the gauge couplings in the low energy limit 
of these theories is invariably a function of the size or shape of the higher di¬ 
mensions. This can be extremely problematic as variation in the gauge couplings 
over cosmological time scales may destroy the successful predictions of primor¬ 
dial nucleosynthesis Recent claims of a time-variation of the fine structure 
constant UDI also motivate the study of theories with dynamically-determined 
(and thus potentially time-varying) couplings. 

In relation to nucleosynthesis, it is only clear what the effect of changing one 
coupling constant at a time has upon the light element abundance. It is possible 
that degeneracies in the effect on nucleosynthesis may arise when more than one 
gauge coupling changes at once: an overall change in all couplings might be 
acceptable at a level much greater than that permitted for any one on its own. In 
this situation it is not possible accurately to constrain models such as j 3 ] where 
the gauge couplings oscillate with a fractional change of the order of 10 % in the 
matter dominated era. Such a large fractional change would not be acceptable in, 
for example, the electromagnetic fine structure constant a at nucleosynthesis, if 
this were the only time-dependent coupling. This paper is an attempt to provide 
an additional constraint, which is independent of nucleosynthesis (but which may 
affect calculations of nucleosynthesis), which suffers as little as possible from the 
problem of relating nuclear forces to underlying theory, and which is sensitive to 
a well-defined combination of couplings. 

Deuterium is only produced during nucleosynthesis, as it is too weakly bound 
to survive in the regions of stars where nuclear processes take place. The fact 
that deuterium is still observed today means that variations in the gauge cou¬ 
pling strengths or other fundamental parameters are non-trivially constrained by 
the requirement that the deuteron be stable at all times after nucleosynthesis. 
The fact that the deuteron is so weakly bound also makes it more sensitive to 
variations in the internuclear force. The strong running of 0:3 at low energies 
means that a change in the coupling strength at high energy is manifested in a 
change in the strong coupling scale Aqcd (henceforth denoted by A), by the usual 
dimensional transmutation arguments. Changes in A in turn lead to changes in 
the internuclear force. 

As recently pointed out by Langacker et al. ini, one also expects changes in 
quark masses and in the Higgs v.e.v. v if gauge couplings are unified at some 
scale. Any viable unified theory should accommodate (if not predict) a mecha¬ 
nism for electroweak symmetry-breaking, which may well depend sensitively on 
SUSY-breaking masses, and a mechanism of generating small Yukawa couplings, 
all of which may have a dependence on the unified coupling. Moreover, A is 
sensitive to all coloured particle masses through threshold effects of RG running. 
These effects introduce a large measure of mo del-dependence, since the correct 
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theories of SUSY-breaking, electroweak symmetry-breaking and fermion masses 
are unknown. One may choose for simplicity to set to zero unknown effects in the 
electroweak, SUSY and Yukawa sectors [12], but this runs the risk of neglecting 
terms which are of equal size or larger than the terms kept in the analysis. 

One might also consider “less unihed” models, with more than one dynamically- 
determined fundamental coupling. The heterotic string dilaton S and volume 
moduli T provide a basic example, where gauge couplings and renormalisable 
Yukawa couplings (for canonically-normalised helds) have a universal dependence 
on S, but may be differently affected by changing T. The greater the number 
of independent quantities considered as time-dependent, the less predictive the 
theory becomes and the less meaningful are any constraints. Here we restrict 
ourselves to estimating the dependence of low-energy quantities in a somewhat 
idealised framework with a single dynamical unihed coupling, corresponding to 
the v.e.v. of a dilaton-like held. 

In the hrst part of this paper we calculate the deuteron binding energy by 
considering meson exchange forces, expressing the relevant parameters as a func¬ 
tions of A and the light quark masses rriu, rrid- The result is rather simple: we 
hnd that the deuteron is stable as long as the ratio rriq/A is greater than a certain 
value, where niq = ruu + md- We perform similar calculations for the dineutron 
and diproton systems in the same isospin multiplet and investigate the criterion 
for their stability. Then we relate A and the quark masses to the QCD coupling 
strength and running masses at high energy using renormalisation group (RG) 
evolution. We take two cases, supersymmetric models with RG running up to the 
GUT scale (similar results will be obtained in the case of power-law unihcation in 
large extra dimensions), and nonsupersymmetric low-scale models with RG run¬ 
ning up to a scale of a few TeV. The main result of this section is the exponential 
dependence of A on the perturbative strong coupling as at high energy. We also 
hnd how the Higgs v.e.v. and SUSY-breaking masses inhuence the low-energy 
parameters. 

Finally we consider how the bounds deduced from the two-nucleon system 
apply to various types of high-energy model, in which as and the quark masses 
depend on model parameters (in particular the sizes of extra dimensions) which 
may be time-dependent. Thus, bounds on the possible cosmological evolution 
of such models since nucleosynthesis can be obtained. The constraints from the 
dinucleon system will in general apply to a different combination of theory param¬ 
eters from those arising from nucleosynthesis — taking the two together bounds 
the variation of fundamental parameters in two directions. (In addition, there 
are many other observational bounds applying at various much later epochs, dis¬ 
cussed for example in ^2), including some from direct laboratory measurement.) 
We also point out for the hrst time and estimate the effect of changing the 
deuteron binding energy on the process of helium formation at nucleosynthesis, 
which may give rise to stronger bounds. 
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1.1 Relation to recent work 


The relation between a time-varying fine structure constant and other observables 
in particle physics was also treated in [iniiaiini. in m the resulting variation in 
the vacuum energy Vq was estimated from general principles of QFT and found to 
be enormously larger than the cosmological bounds; hence, the authors concluded 
that a large number of implausibly accurate fine-tunings would be necessary for a 
time variation of the size that has been claimed to be consistent with field theory 
and cosmology. Such an argument is rather weak since it assumes, crucially, that 
the cosmological constant problem is “solved” at the present time by cancellations 
of different field theory contributions, all of which are many orders of magnitude 
larger than the measured value of Vq. This can hardly be a solid starting-point 
from which to set theoretical limits. The alternative considered in ffH was a 
self-tuning mechanism which protects the four-dimensional spacetime were we 
live by dynamically “absorbing” the vacuum energy into the curvature of one or 
more extra dimensions. Such a mechanism probably disallows inflation, but it 
is by no means clear that it also rules out spontaneous symmetry-breaking as 
claimed (the argument being that anything that “cancels” the vacuum energy 
also removes the source for a scalar to roll to the minimum of its potential). The 
existence of a dynamical time-scale for the self-tuning, the role of thermal effects 
in creating an effective potential, and the possibility that the extra-dimensional 
model may break four-dimensional Lorentz invariance are possible ways out of 
the D = 4 field theory argument that V (0) must vanish at all times (see ^3] for 
related discussions). 

We prefer to take the majority point of view on the cosmological constant, 
i.e. that it is our greatest theoretical area of ignorance and that no credible 
way to explain its smallness currently exists, therefore very little can be deduced 
from it, and certainly nothing related to quantum effects in the theory (which 
are the main difficulty in accommodating varying alpha). If we want to retain 
the semiclassical picture of matter coupled to gravity, the only sensible interpre¬ 
tation is a very light scalar, which interestingly would have about the same mass 
as quintessence; due to the coupling to electromagnetism, such a scalar would 
mediate composition-dependent forces which might have experimental signatures 
ra- But the vanishingly small scalar mass (even if experimentally confirmed) 
would remain a mystery, in the absence of an underlying theory which would ex¬ 
plain why spacetime curvature was apparently so insensitive to a cosmologically- 
evolving field theory. 

Calmet and Fritszch ^2] calculate some of the consequences for low-energy 
physics of changing a, within a GUT-like theory which constrains the SM gauge 
couplings to be equal at a particular energy scale (or at least to satisfy some fixed 
relation). They consider exclusively the effects on the strong interactions, with 
the unstated assumption that the mechanisms of electroweak symmetry-breaking, 
supersymmetry-breaking and fermion mass generation are held constant despite 
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varying the unified coupling, hence that the quark and lepton masses, as well 
as W and Z masses, remain unchanged. Our calculation of the dependence of 
A and Mjv on 0 ^ 3 (p > rrit) is essentially identical to theirs, except for including 
the full dependence on thresholds. As noted above, such a study can only be a 
first step since the electroweak mass scales may also depend sensitively on the 
unified coupling (for example in models of hidden sector SUSY-breaking, radiative 
electroweak symmetry-breaking and anomalous U(l) flavour models). 

Unfortunately, since we do not know the correct theories of electroweak and 
supersymmetry-breaking, let alone that of fermion masses, correlations between 
the time-variation of different low-energy observables (such as a and the proton 
magnetic moment /ip or the mass ratio irtp/me) are necessarily model-dependent, 
unless one can find a combination insensitive to, say, the electroweak sector. 
We take the converse approach and quote in a hopefully less mo del-dependent 
way the bounds on high-energy parameters deriving from the low-energy system 
that we are studying. In particular one cannot claim yet (as in M) that the 
current data on the time variation of a and other quantities are inconsistent 
with unification; but given an observational bound one can derive bounds on 
time-varying fundamental parameters, in whatever model one is interested in. 

With precise measurements of at least two different quantities at any partic¬ 
ular epoch, of which at least one shows a nonzero variation, one can rule out 
classes of unified theory that predict the wrong relation between (variations in 
the) different quantities. Time-varying couplings in principle are a new way of 
doing phenomenology, which allow us to test relations between the derivatives of 
different quantities rather than their static values. This paper presents an addi¬ 
tional bound which applies at all times after nucleosynthesis, thus for models in 
which time variation was more rapid in the early Universe (such as ^) it is likely 
quite restrictive. 

Langacker et al. ^2] have a somewhat similar approach to the “high-energy” 
aspect of the problem, parameterising the effects of variations in the unified 
coupling strength on the electroweak and Yukawa sectors by unknown, model- 
dependent constants of proportionality which depend on the particular model 
used. By looking at other precision astronomical measurements besides mil. they 
find that these constants must satisfy a rather precise relation, for the model to 
be consistent with observation. 


2 Nuclear forces and the stability of di-nucleons 

2.1 Chiral symmetry breaking and the pion 

Previous calculations taking into account the effect of the time variation of the 
strong interaction on low energy nuclear phenomena were performed in the chiral 
limit, i.e. the limit of massless quarks and pions, where the only energy scale of 


5 


the system is A. The approach of was to assume that certain dimensionful 
static quantities (for instance vector meson masses and nuclear binding energies) 
were directly proportional to A; in this way it is possible to place constraints on 
the variation of over time. 

Explicit breaking of chiral symmetry in QCD is achieved by the addition of 
non-zero quark masses to the Lagrangian. The spontaneous breaking of chiral 
symmetry occurs when the quark condensate (0|gg|0) develops a non-zero value 
and dynamically lines up with the quark masses in the internal SU(nj) flavour 
space The pion is the Goldstone boson associated with this spontaneous 

symmetry breaking: it is not quite massless, since the chiral symmetry is explicitly 
broken by the quark masses. The pion is however much less massive than any 
other strongly-interacting particle, which means that the internuclear pion force is 
relatively long range and important in the analysis of the loosely bound deuteron 
and unbound di-neutron and di-proton. 

Models of nuclear structure utilising just the exchange of the a;, a and p 
mesons jSS] are quite successful. This is partly because the isospin dependence 
of the pion interaction is such that the average effect across an isospin symmetric 
nucleus consisting of many protons and neutrons cancels. However, no such 
cancellation occurs in two-nucleon systems, where the internuclear pion force 
is hnite and well-defined. The large separation of the proton and neutron in 
the deuteron (and in the putative di-neutron and di-proton systems, if coupling 
strengths were to change so that they were weakly bound), relative to the typical 
range of internuclear forces, increases the importance of the long range pion 
force in the binding of the system. Consequently, the pion force is the dominant 
contribution to the binding of two-nucleon systems, and the contributions of a;, 
(j and p meson exchange to the binding energy can be taken to be of secondary 
importance. 

The pion mass is given by the Gell-Mann-Oakes-Renner relation [21] 

= i.'rnu + md){0\qq\0) (1) 

The non-zero pion mass leads to a hnite divergence of the total axial current, 
which compensates for the partially conserved axial current of the weak interac¬ 
tion Hamiltonian. This compensation leads to the Goldberger-Treiman relation 
for the pion-nucleon coupling Qt, |21j 

‘^MnQa 

9n = —7 - ( 2 ) 

In 

where Mjv is the nucleon mass. Although this is an approximate relation, one 
would expect it to remain valid as the parameters fT, and change, so long as 
the variation is much less than that necessary to restore the chiral symmetry of 
the vacuum at zero temperature. (If this were not the case, the problems created 
for nuclear physics in the early Universe would be much greater than those we 
discuss here, indeed such a scenario would be immediately ruled out). 
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Since the nucleon mass Mjq originates mainly from confinement of the quark 
colour charges (rather than current quark masses), Mjv simply scales in direct 
proportion to A. Combining this with Eqs. o and 0 we hnd the relation 

glm^ + md)^/2(0|gg|0)^/2 gl{mu + 

Mn ~ 2gAMl ^ gAh? ^ ’ 

Note that the product mq(0|gg|0) is invariant under change of renormalisation 
scale in QCD (see e.g. j2n!) so thus far we avoid the problem of having scale- 
dependent quantities in a context where mass scales may be time-dependent (see 
section EH). 

To proceed further, we need the relationship between (0|gg|0) and A, obtained 
by looking at the effective potential for the chiral symmetry-breaking, which 
shows that the energy scale associated with chiral symmetry restoration varies as 
the value of the condensate: 

(0|gg|0)5~A, (4) 

a result also predicted by sum rule calculations of the nucleon mass |2H1. As it 

stands this is not a well-dehned relation, since the LHS is renormalisation-scale 
dependent. To remedy this, we use the formalism of “invariant quark masses” 
described e.g. in to one loop the equation for the running quark mass mi{g) 
is solved by 

mi(p) = mi(-/3ia3(/i)/7r)“'^A/3i (- 5 ) 

where rhi is a RG invariant quantity, j3i is the one-loop beta-function coefficient 
and 7 i the one-loop anomalous dimension of quark masses, such that for three 
flavours — 7 i//di = —2/(—9/2) = 4/9. Thus for scales below rric we have = 

up to a constant universal factor which will drop out of fractional 
changes A(lnm) = Am/m. Then the bilinear order parameter must have the 
one-loop behaviour 

(0|gg|0)(/i) = (RG invariant) x 0:3 (p)“"^/® oc A^a 3 (p)“A 9 (^ 5 ) 


which gives the correct dependence on both the RG scale and (a possibly time- 
varying) A. We can now write 


Mm 


2 I mu+ md 
const. X g^ ( -- 


(7) 


in which all RG dependence cancels neatly We will see in the following sections 
that the value of the ratio (m^/A)^/^, which we will denote as c. 


c 


2 _ 


mg 

A 


^Such a redefinition can be performed to any desired loop order. 


( 8 ) 
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is the parameter controlling the effect of varying coupling strengths on nuclear 
physics phenomena dominated by pion exchange. Later we shall see how c can be 
related to coupling strengths at high energy in some examples of unihed models. 

Many different estimates of the dependence on nuclear binding energies on the 
one dimensionless parameter exist; despite the fact that the deuteron is the 
simplest nuclear system, a level of understanding sufficient to estimate reliably 
from hrst principles is lacking. A systematic approach based on expanding in the 
light quark masses has been proposed cni, but even in this framework the lack of 
control of four-nucleon operators introduces uncertainties in the estimation of Bd 
(highlighted in jT^j, which appeared after the hrst version of this paper). Since 
we do not claim to be doing precision calculations of Bd, and would be satished 
with an estimate of its dependence on which was inaccurate by a factor of 
a few (remember that we are placing cosmological bounds on 6a), this is not a 
major concern. However, Beane and Savage claimed in HZI that the dependence 
on could be mostly or entirely erased for some values of coefficients of these 
four-nucleon operators, drawing the conclusion that bounds on varying couplings 
from nuclear physics were considerably weakened. Without further analysis we 
cannot tell to what extent this particular choice of coefficients is hne-tuned, but in 
the absence of an underlying reason it appears very unlikely that the ultraviolet 
effects (from the point of view of the chiral expansion) giving rise to the quark 
mass-dependence of four-nucleon operators would conspire with the low-energy 
effects of pion exchange in such a way. Since all previous estimates gave a rather 
steep dependence of Bd on c^, such a conspiracy would be a very unlucky (or 
lucky?) coincidence. 


2.2 The di-neutron 


In this section we attempt to calculate how large a change in the parameter c 
would allow two neutrons to form a stable bound state. The total spin of the di¬ 
neutron ground state would be zero and the neutrons share the same orientation 
in isospin space so we use the potential for the S' = 0, / = 1 state 1221 

f2 p-rriTrr 

^(r)(5=o,/=i) = -f ^- (9) 

An r 

where the dimensionless coupling is given by 


r 9i (y 

An An \2Mjsf) 


7.95 X 10“2 


( 10 ) 


where the numerical value is derived from present-day measurements. We assume 
a trial wavefunction of the form 


^(r) = (,-^lm.T^-hvn^T 


( 11 ) 
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Figure 1: Total ground state scaled energy of the di-neutron with respect to 
the parameters b and v introduced in the text. The light region corresponds to 
positive energy (unbound) and the dark region negative energy (bound). 


In this equation -^(r) = rT(r), where T is the radial part of the wavefunction, 
such that 'ip‘^{r)dr is the probability of finding the nucleon separation to be be¬ 
tween r and r -|- dr. At large r where the Yukawa potential is negligible, the 
Schrodinger equation for the relative motion of the two nucleons gives 


1 d^-^(r) 
Mn dr'^ 




( 12 ) 


where Eb is the (negative) binding energy of the state and the wavefunction 
m takes the form 


^/>(r) = 


\/—EbMjsi 


(13) 


The trial wavefunction m is a suitable choice since at small r it is independent 
of the binding energy whereas at high r it is completely determined by the value 
Eb- By applying the Hamiltonian to the wavefunction (fTT|) we hnd the energy of 
the system ins 


H(r)'ll)(r)dr 

rmrWdr 

x/6/^m^/Co(2v/2T46) b^/^mlIC2{AVb) 
47r/Ci(4v^) ^ 2MN]Ci{4Vb) 
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where the fCn are modihed Bessel functions of the second kind. We dehne a 
dimensionless parameter v which is directly proportional to c dehned above 


P 2Mn 
An m-it 


oc c. 


(15) 


Thus, we can express the scaled ground state energy, in units of 6m^/2M7v/Ci(4v^), 
in the simplihed form 

-^'scaled 6/Co ( 2 P 2 + Ab^ - n/Ci(4). (16) 

The dependence of T^scaied on b and v is illustrated in Fig. [T] In order for a bound 
state to exist, we require v > 2.8, or 


2 An 47r M^r 

We denote by vq the currently-observed value of 
1.08. Then the di-neutron stability criterion is 

Esinn) < 0 ^ — = -^ > 2.6 

Vo Co 


(17) 

numerically equal to 

(18) 


thus if c = prhq/K increases by a factor of 2.6 the di-neutron will become bound. 


2.3 The di-proton 

The fact that nuclear forces are independent of electromagnetic charge is illus¬ 
trated by the result that the difference between the binding energies of and 
He^ is well explained by the energy of Coulomb repulsion between the two pro¬ 
tons in He^ EH- One would expect the ground state of the di-proton to have the 
same nuclear quantum numbers as the ground state of the di-neutron. We can 
therefore see from the previous section that reducing the Coulomb repulsion to 
zero will not be enough to bind the di-proton: only a large variation in the strong 
force (more precisely, in the ratio rhq/A) would achieve this. 

One would expect the size of a marginally bound di-neutron or di-proton to 
be of the same order as the effective range of the potential responsible for the 
binding. For a Yukawa potential like the one in Eq. o that we are using, the 
effective range is of order of the pion mass 


Ceff ~ — (19) 

m.^ 

and the Coulomb repulsion energy at that range is approximately 1 MeV. Thus 
the parameters necessary to obtain a bound di-neutron with a binding energy of 
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1 MeV will be similar to those required to bind the di-proton at threshold. We 
assume here that the effect of varying a on the electromagnetic repulsion will be 
dwarfed by the corresponding change in Trig/A. This assumption is shown to be 
justified consistent in the next section. Equation (|T2jl and Fig. [T] show that this 
corresponds to the condition on the parameter v of 

Eb(pp) < 0 — = — >3.2 (20) 

Vo Co 

so if c increases by a factor of 3.2 the di-proton will become bound. 

There also exists a hypothetical proton-neutron pure L = 0 state which one 
would expect tro be stabilised for similar values of c as the di-neutron. 


2.4 The deuteron 


The deuteron has total isospin / is zero, so the spins of the two nucleons must be 
parallel by the Pauli exclusion principle. The L = 0 wavefunction is unbound in 
the same way as the di-neutron and the di-proton, however the parallel nucleon 
spins mean that the L = 2 wavefunction is bound by the tensor force and the 
ground state forms an admixture of these two wavefunctions. The exact form 
of the wavefunctions is given by the solution of a coupled pair of 2nd. order 
differential equations 122 ]. To solve these equations explicitly for different values 
of the parameters /, tUtt and Mn is highly nontrivial, so we will resort to a simple 
square well model of the kind used in |27|. 

The analysis is then simply the bound state condition for a finite square 
well potential of depth Vsq and width a, found (assuming \Eb\ |Wg|, which 
is the case in the situations considered) by matching a sinusoidal solution for 
the wavefunction inside the well to an exponentially decaying mode outside. We 
assume a = m~^ then the binding energy Eb is then given by the solution of 


cot 




( 21 ) 


Then, using the observed values of the pion mass and the binding energy Eb = 
—2.226 MeV we find that the depth of the well is Vgq = —66.15 MeV. Since the 
prefactor in the internuclear pion potential P/An is dimensionless we need to 
see how the quantities in this simple model scale with the parameter c that we 
have been using to parameterise the variation in the underlying gauge couplings. 
The corresponding dimensionful ‘depth’ Vy of a Yukawa potential is related to 
equation O by 


Vir) = V, 


-rriT^r 


m^r 


YjL 


f 


2 


4n 

P 


4n 


r 


( 22 ) 
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so we can expect Vsq, the depth of our square well, to be directly proportional to 
pm.„. Making the substitution /3 = Es/Vsg equation (pT|l becomes 


cot 


MjyVsq 

ml 


(/9-1)) 



so the relevant combination of parameters is MjqVgq/ml 


M^Vsq 

ml 




Mn 

m^ 


oc c. 


(23) 


(24) 


Using Eq. (1231) we can see that the value of c where the deuteron becomes unstable 

is 


^unstable 

Co 4 MN\Vsq 


(25) 


so a 13% reduction in c will give rise to an unbound deuteron. 


3 Running of gauge couplings and unification 

3.1 GUT-like models 


The LEP precision measurements of gauge couplings suggest that the three stan¬ 
dard model gauge groups become unihed at some high energy Mjj ^ 2x 10^® GeV 
j2j. Thus any variation in the unihed gauge coupling ajj{Mu) leads to calculable 
changes in the low energy values of the gauge couplings. The dependence of low 
energy couplings on au is very similar in the case of large extra dimensions jH] 
or non-factorisable geometries jHj where the unihcation scale is much lower. This 
is because extra contributions to the renormalisation group equations from the 
massive Kaluza-Klein modes of the theories change the running in such a way 
that the effective energy scale and coupling strength of unihcation stays the same 
from the point of view of the low energy 4D theory (see, however, EO). In the 
following analysis we start with the simplest scenario of a “grand desert” between 
the weak and GUT scales with M = 1 supersymmetry. However, no matter what 
the matter content or other details, results of a similar order to those presented 
here seem unavoidable for any theory with gauge unihcation. 

The running of the gauge couplings Oj, i = 1, 2, 3, with the renormalisation 
scale p is given at one loop by the expression 


dai 

dt 



t = ln(/^//^o) 


( 26 ) 


where po is an arbitrary, constant reference scale and the bi depend upon the 
gauge groups and matter representations transforming under each group. In 
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SU(5) and SO(IO) SUSY GUT’s and a large class of string models the gange 
conplings are nnified at some high energy scale Mf/: 


— a2{Mu) — (5/3)ai(Mt/) — au{Mu) (27) 


where the normalisation factor 5/3 derives from the fact that the U(l) hyper¬ 
charge gange gronp mnst be embedded consistently in the nnihed gronp or string 
model. Other normalisations are possible, bnt are often inconsistent with nnih- 
cation [Hj. At energies above the snperpartner masses, the 63 and 62 of Eq. (jSED 
are given by [7] 

&3 = 9 - 2nG, 62 = 6 - 2nG - ^ (28) 


where Uff is the nnmber of Higgs donblets and no is the nnmber of fermion 
generations. Integration of Eq. (Ei above the scale of snpersymmetry-breaking 
yields 


Mu 

dvr 

( 1 

1 V 

M "’'P 

_6 -h uh 

1«2(M) 

a,{M)J\ 


(29) 


Using the standard relation ai cos^ 9w = ^2 sin^ 9w = « and snbstitnting the 
observed valnes of siv?9w-i (^z(Mz) and a(Mz) |H] in this eqnation, we hnd a 
nnihcation energy Mu 10^® GeV and a nnihed conpling of ~ 25. 

In order to hnd the ehect on the strong interactions of changing au we will 
nse the 63 of Eq. (El to rnn from Mu down to low energy for diherent valnes 
of au- Then we assnme that all the snperpartners of the standard model particles 
have the same mass m, which dehnes the scale below which snpersymmetry is 
broken. At energies less than rfi the rnnning of 03 proceeds as in normal QGD, 
i.e. 


1 _ 1 _ 33 - 2nf / Mio„ \ 

a3(A'4ijh) a3(A/|„„) 6?r 


(30) 


where is the nnmber of qnark havours of mass rriq < mhigh- We nse the observed 
valnes of the top, bottom and charm qnark masses, dehned at the self-consistent 
renormalisation scale mq(mq), in this part of the rnnning. 

If one neglects qnark masses, the Lagrangian of QGD contains no dimen- 
sionfnl parameters: however, on qnantisation, loop corrections give rise to the 
renormalisation gronp invariant^ strong interaction scale A: 


'' ^'^“'’((33-2r!/)a3(A4) 

where M < rh. This scale sets the characteristic energy of the particles of the 
low energy effective theory. Even in the presence of massive qnarks, this qnantity 

^Invariant unless one crosses the mass threshold of a quark or snperpartner species, which 
changes the running and the apparent value of A. Here, we calculate A under the assumption 
that it lies at the energy scale rug < A < rric- 
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Figure 2: Variation in A (strictly, A/Mu) vs. variation in the gauge unified gauge 
coupling. 


remains important in understanding low energy phenomena, as we saw in the 
calculation of internuclear forces. 

For constant au, the running of below the superpartner thresholds and 
the value of A/Mu depend strongly on the common superpartner mass m. But 
if fh/Mu is held fixed and au varied, the fractional change in A/Mu for a given 
fractional change in au has little dependence on the value of m. For any constant 
value of fh/Mu correponding to superpartner masses between 100 GeV and 2 TeV, 
we find the same ratio of fractional changes {Mu/A)A{A/Mu)/{Aau/au), to well 
within the accuracy of our results. Later, we give the explicit formula for A in 
the case fh/Mu is also changing. 

In order to obtain a^ at low energies and A/Mu as a function of au{Mu) 
we include the masses of the top, bottom and charm quarks jS] and change 
the running accordingly, with results illustrated in Fig. |21 To first order in the 
fractional change Aau/au, we should be able verify the results for A (In A) in 
nacsi. Later, we will consider the effect of varying the threshold masses, but 
for the moment we consider them as fixed. Using the result quoted later as Eq. 1221 
we find 

Gtt 

A(lnA/M;7) = —a/j^A{\nau) ~ 17A(lnat/) (32) 

consistent with the result of both papers. (Note that any determination of the 
absolute value of A/Mu will depend significantly on the renormalisation pre¬ 
scription, so without a detailed treatment, including also higher loop effects, our 
estimate of the prefactor is likely rather imprecise.) 
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Figure 3: Variation in A/M* vs. variation in as at a scale of 1 TeV. 

3.2 Non-GUT (brane world-type) models 

One can perform a similar analysis in the case of theories where the gauge cou¬ 
plings at the fundamental scale do not satisfy a GUT relation, for example inter¬ 
secting brane models (see e.g. uniiisi) where different gauge groups propagate 
along different world-volumes (|11]). In this case we just hnd the relation be¬ 
tween the SU(3) coupling just below the fundamental scale M* and A, by the 
same method. For a high-scale model with softly-broken SUSY the calculation 
is identical, replacing au by ~ M*) (where M* is the scale below which 
d = A effective theory applies for gauge interactions). In the case of a low fun¬ 
damental scale of order 1-10 TeV and no supersymmetry we simply run from 
M* ~ TeV down through the heavy quark thresholds, with results as illustrated 
in Fig. El The main lesson from this exercise is the sensitivity of A/M to the 
strong coupling at high energies M S> Mjy: the slope of both graphs around 
A = 0 is 10 to 20. In our conclusions, for simplicity and to avoid making highly 
model-dependent statements, we will pretend that the analysis hnishes at this 
point, and express the results in terms of a fractional change in au under the (in 
general unrealistic) assumption of constant quark and superpartner masses. 


4 Thresholds and RG running 

In both cases (GUT-like and non-GUT) the running of as is complicated by 
the dependence on heavy quark and coloured superpartner (gaugino and squark) 
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thresholds. In the case with superpartners the one-loop result with tree-level 
matching at thresholds is 


A 


]^g-67r/27a3(M) 


/ 3 _ 3 \ 2/27 

/ mcTTibmtm'lmq ' 




M9 


(33) 


where m\ is the gaugino mass, rriq is the geometric average squark mass and M 
is an arbitrary scale above the superpartner masses; for the SM we find simply 

(34) 

where M' > rrit. The quark and superpartner masses are the “decoupling masses” 
mq{^ = rrtq) etc. in a convenient renormalisation scheme. 

The possibility of varying threshold masses complicates the simple one-loop 
formula 

1 Aa(/i) 1 Aa{fi') 

a{fi) Q;(/i) «(/!') 

derived in [T2] for a hxed ratio if there is, say, a quark with /i < < /i'; 

we hnd 

1 Aa;(/i) 1 AQ((/i') ^ ^ j mq^{mq^) 

a{fi) a{fi) Q((/i') Q;(/i') 27r mu' 

where simply denotes the contribution of qi to the beta-function 

coefficient above its threshold. 

One expects higher loop effects to change these relations, generically by the 
introduction of a power-law dependence of A on (etc.!) 


4.1 RG scale ambiguities 

What is meant, in the context of a time-varying theory, by quark masses “at Mz" ? 
Since in all likelihood the Higgs v.e.v. and the SU(2) coupling will also be time- 
dependent in a unified theory, such an expression becomes ambiguous through 
the change in Mz- It is commonly noted that dimensionful quantities cannot 
meaningfully be said to be time-dependent, since one has no way to guarantee that 
units of measurement do not also change; the presence of RG scale dependence 
threatens to introduce dimensionful quantities through the back door unless we 
take care in our notation. 

One convenient approach is to dehne a reference scale which can always be 
set to a constant without loss of generality; in discussing unified theories the 
obvious choice is Mu (or we may choose the Planck scale, or the string scale, 
or the “quantum gravity” scale M* according to the variety of theory we are 
considering). Then we would write in all explicitness 

mq{n/Mu = Mz/Mu = f{oiu)) 
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remembering that Mz will depend in the unihed coupling which itself is time- 
dependent; alternatively one could set a hxed ratio equal to the current 

value {Mz/Mu)\o ~ 4 x 10“^^ to avoid the complication of calculating Mz/Mu 
as a function of au. Once we realise that the RG scale can thus be converted to 
a dimensionless parameter, the ambiguity disappears. 

Now we tackle the case of time-varying decoupling masses for heavy quarks 
in more detail; let us take the bottom quark for dehniteness. The running mass 
is given at one loop by 

rribifi) = rhi,{a2,{mu)y^^^^ 
where the bottom mass runs as 

—mb = -ruftl—as H-) 


and 7 i = 4 in QCD (note that hi = 11 — 2?7,//3 in our notation). The decoupling 
mass mf satishes 


1 ,, 6i , mi 
an^iM) -F — In—^ 
^ ^ ^ 27r M 



-61/71 


where M is a constant reference mass above mb but below mt and we have sub¬ 
stituted the one-loop solution for as on the LHS. Differentiating and rearranging 
to isolate A In we hnd 


(1 + Aln(m^/M) 


A In(m6/M) -f- 


7i azimf) 
hi aii{M) 


Alnas(M). 


(36) 


As remarked before, the variation in ifib/M is given just by the change in the 
Higgs v.e.v. and in the Yukawa coupling yb at high scale (since we are at the 
moment only considering the QCD contribution to the running of quark masses, 
which will be the dominant effect at low energies). 

Note that the second term in brackets on the LHS of Eq. is formally of 
higher order than the hrst, therefore we may consistently discard it for the heavy 
quarks for which as is perturbative. The second term on the RHS is also a factor 
of as down compared to the analogous expression for the variation of the QCD 
scale A (see Eq. OHl), but we cannot tell without knowing the dependence of mb 
on au which term on the RHS may dominate. If the second term on the RHS 
does dominate, then the changing strength of QCD is the main cause of mf/M 
varying, in which case quark threshold effects will be a rather small correction 
to the change in A and can be neglected. If the hrst term dominates, the main 
effect is through the variation of quark masses at high energy, in which case we 
can take Aln(m^/M) ~ Aln(mfe/M); this is the relevant formula if one expects 
quark mass thresholds to play a signihcant role. 

We also note, on the subject of fermion masses, that since the hne structure 
constant is dehned at the scale m^, a variation in m^/Mu can induce a change in 
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a even if au is unchanged; this was pointed out by Wetterich j31j in a preprint 
appearing after the hrst version of this paper. However due to the very slow 
running of a{fi) one would require a rather large change in rrie/Mu, or in any 
other charged particle threshold, to reproduce the observed variation in alpha. 
This restriction, taken together with the bound on variation of Mp/irte ISSI, imply 
that a sizable contribution to varying alpha from this kind of threshold effect is 
unlikely 


4.2 Dependence of c on high-energy parameters 


Now we attempt to express (? = {mg/A) as a function of some more fundamental 
parameters. Substituting for mg = yv/\/2, where y is the average light quark 
Yukawa coupling, modihed by the appropriate function of 03 (/i), and using the 
expressions derived above for A we hnd 


m = const. X 


in the SM and 


c = const. X 


M{yc{yc)yt{yt)yt{yt)!ilcyvV^rniy'‘/^^ 


(36) 


(37) 


in the MSSM, where S/ 3 {cy) = sin(cos)/3. From now on we assume that the 
change in tan/3 is small and lump the light quark Yukawas together as yq. Since 
the running of the up, down and charm Yukawas as they appear in this expres¬ 
sion (as-corrected for the light quarks) over the range up to Mu ~ 10^® GeV is 
small, we set them equal to their value at unihcation yq{M^). In fact, to hrst 
approximation we neglect the nonlinear running of the b and t Yukawas also: in 
a careful calculation one would use the semi-analytic solution for the top Yukawa 
and include the effects of the top feeding into the other masses, but the effects 
that we are neglecting by effectively taking yt{y) = yt{M^) for y, > mt are likely 
to be sub-leading compared to, say, a large fractional change in the Higgs v.e.v. 
or in SUSY-breaking masses. 

Thus we reach the approximate expressions (which still likely contain the 
leading dependence on au) 

7 

A(lnc2) = p^A(lna 3 (M[ 7 )) +-A(lnu/M[;) + A(ln//g[;) 

2 2 

- —A{\nycu + liaybu + ^^Vtu) - -A{\Yimx/Mu + \nfn^/Mu){38) 
for MSSM unihcation and 

27r 7 

A(lnc^) = p^A(lna 3 (Md) +-A(lnu/Md + A(ln//gd 

2 

- — A(ln//c* -F Inyb* ln//f*) (39) 
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for the SM running in low-scale models. 


5 Finding bounds on fundamental models 

If we assume that the supersymmetric model we are referring to is associated 
with the low energy effective theory of some heterotic superstring or M-theory 
compactification the parameters au and Mjj are related to the size of the six 
dimensional manifold upon which the theory is compactihed.^ In this situation 
au oc V~^ oc {MuY where V is the volume of the compactihcation manifold [Hj. 
Because of this we might expect that as one changes the size of the extra dimen¬ 
sions, both the GUT unihcation energy scale and coupling might vary together. 
However the effect of the changing energy scale has a negligible effect on the 
results because of both the 6 th power in the relation between the energy scale 
and the coupling, and the slow running of the gauge couplings at high energies. 

The supersymmetry-breaking masses entering into the expression (P?7|) depend 
on the mechanism of SUSY-breaking, which can either take place through non- 
perturbative effects in a hidden sector (as is appropriate in models with high 
fundamental scale) or perturbatively by a Scherk-Schwarz mechanism in mod¬ 
els with large (TeV“^-size) extra dimensions Taking the scenario 

of gaugino condensation in a hidden sector and string moduli mediation as a 
benchmark, the functional dependence on the unihed coupling is 

fh oc 7713/2 OC A^Mp'^ (40) 

where Mp is the reduced Planck mass 1 /a/SttGat and the hidden sector conhne- 
ment scale is ~ Thus the fractional change expected is 

related as 

Gvr 

AiXnm/ Mu) — 2A log(M[//Mp) — --A(lnQ!; 7 ) (41) 

hnotu 

where the hidden sector gauge coupling is taken also to unify at Mu- Using 
the relation Mu ~ Mg oc guMp from heterotic string theory, and imposing the 
correct magnitude of ~ fh/Mp ~ 10 “^® we hnd 

A(ln?7i/Mp) ~ (2-|-34.5)A(lna;7) (42) 

where the factor 34.5 ~ 15 In 10 is universal to models of SUSY-breaking by strong 
coupling effects in a hidden sector mediated by non-renormalisable operators (as 
also found in [TH]b 

^Although the low energy limit of M-theory is 11 dimensional, the gauge degrees of freedom, 
at least in the strong coupling limit of the heterotic theory, are confined to a 10 dimensional 
brane. The appropriate volume to consider is then that of the 6d Calabi-Yau manifold, as in 
heterotic string theory. 

"‘if the fundamental scale is intermediate (10‘‘-10‘'^ GeV) then SUSY-breaking may occur 
directly at tree level in a non-supersymmetric hidden sector and be gravitationally-mediated, 
a possibility recently realized in certain D-brane models m 
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The quark Yukawa couplings gain a universal factor of {S + S) ~ gu 
from the normalisation of the superpotential in going from SUGRA to softly- 
broken global SUSY, thus in particular ytu will vary proportional to ajj which 
will be of importance for electroweak breaking. The light quark Yukawas are 
extremely model-dependent, but estimates can be made in some classes of uni¬ 
fied models. The most common mechanism employed to generate small Yukawa 
couplings is the Froggatt-Nielsen picture in which effective couplings arise from 
non-renormalisable operators when scalars charged under a U(l) group get v.e.v.’s 
{X). In fact this picture finds a natural embedding in the heterotic string where 
the group is now “pseudo-anomalous” and is broken near the string scale due to 
a nonzero Fayet-Iliopoulos term at one loop destabilising the symmetric {X = 0) 
vacuum. The fact that this is a one-loop effect tells us the dependence of the 
v.e.v. on the unified coupling as {X)/Mp oc gp (see e.g. ISHI), thus the small 
Yukawas are generated schematically as 

y,u ~ {Xf^ oc g^^ (43) 

where the quark g* has a U(l) charge Qi of opposite sign to Qx- Thus with 
knowledge of the charges, which also follow phenomenologically if we know the 
constant of proportionality of {X)/Mp = kgp and demand the correct values of 
yiu, one easily finds the variation in ym as a function of A(lna[/). 

Finally, as is well known the standard scenario of radiative electroweak sym¬ 
metry-breaking is sensitive to the soft masses rh and the top Yukawa ytu (for a 
review, see EH). While the full RG equations for the top Yukawa and the rele¬ 
vant soft masses are complicated, the leading dependence on the input quantities 
was already considered in |1H] in which the following estimate was made for the 
sensitivity of the W mass to yt 

Mw{y<Qo) ^ ytmt\n{y^/Ql) 

^X{\nMw{y)/Mu) ^ A(lnm/Mf;) + A(ln|/,)ij^^^ (44) 

In Qo/At 

where Qq is the scale at which the up-type Higgs mass-squared crosses to negative 
values, given by , where A; is a ratio of SUSY-breaking masses at the scale 

Mu and we choose a fixed renormalisation scale p ~ Mw\o since (analogously to 
the case of heavy quark thresholds) the running of g 2 {g) is slow. Then we impose 
Mu/y — 10^^: one finds from a more careful analysis |1H1 that Qo//^ tends to be 
of order 10, in which case the enhancement factor \a{Mu/g )/ ln(Qo/ g) is of order 
15 (although this also depends strongly on the pattern of Higgs mass parameters 
at unification). 

Hence in addition to the strong SUSY-breaking mass effect previously dis¬ 
cussed, there is a milder dependence of Mw or v through yt oc and the 
overall effect on the Higgs v.e.v. may well be the largest contribution to any 
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low-energy physics that depends on quark masses, even relative to the exponen¬ 
tial dependence of A on the high-energy coupling strength in QCD. The reason 
for this is subtle: while the beta-function coefficients giving the running of the 
hidden sector gauge coupling are assumed constant, thus the hidden sector scale 
“feels” the whole of the hierarchy between Mu and Mw, the presence of 
thresholds at low energy (<C Mu) in QCD means that Aqcd feels mainly the 
running at low energies (< rh) thus the sensitivity to varying the coupling at 
energies above the thresholds is less. If there were also thresholds for charged 
states in the hidden sector at masses Ah ^ Mh Mu, the leading dependence 
of SUSY-breaking on the unihed coupling would be smaller, however such a case 
is atypical for the vector-like matter representations usually considered. 

However, since SUSY-breaking and electroweak breaking are highly model- 
dependent, we will not go through the resulting bounds in detail. Other types of 
model are likely to give completely different results, for example Scherk-Schwarz 
SUSY-breaking, where the electroweak and SUSY-breaking mass scales vary in¬ 
versely with the radius of an extra dimension, as do the 4-d gauge couplings 
(more precisely, Oj oc R~^) and so the fractional variation in v would likely be of 
the same order as that in au (rather than maybe 50 times larger). 


6 Discussion and Constraints 

In drawing our conclusions about the bounds on the variation of au, we will 
simply for the sake of an easy comparison set to zero all variation in quark masses 
and other thresholds, from no matter what source, leaving the effect of on A as 
the only varying quantity. This is not intended as a realistic treatment of unihed 
theories but only as an illustration of one effect, which can be systematically 
extended to whatever unihed theory by going back to the bounds on c = (m^A)^/^ 
and substituting for the expressions we derived in Sections |3] and El With this 
caveat we proceed to discuss the bounds from the di-nucleon systems. 

6.1 The deuteron 

A 69% increase in c~^ = A/(m„ -|- rrid) will destabilise the deuteron. On the 
assumption of hxed quark masses. Fig. |21 tells us that this corresponds to a 3% 
increase in the gauge coupling at the scale of unihcation for a conventional SUSY 
GUT. This creates problems for cosmological models where the gauge couplings 
vary signihcantly over cosmological time scales jl], since no such variation can 
have occurred at any time since nucleosynthesis. 

For a model with a TeV scale GUT, 0:3 cannot have increased more than 
8 % at any time since nucleosynthesis, since this would also lead to the deuteron 
becoming unbound. 

Recent observations suggest a negative variation in the electromagnetic hne 
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structure constant at high redshift HU, which would, if one assumes high scale 
gauge unification and consequently a corresponding negative change in a^, in¬ 
crease the binding energy of the deuteron. Unless the field responsible for the 
value of the gauge coupling is frozen by the Hubble expansion HE] one might 
expect this variation to be much larger in the early universe (see e.g. [HHI)- This 
might result in the deuteron having been more resilient to photo-dissociation by 
the decay products of massive relic particles [211 • This possibility would be ben¬ 
eficial to models where the decay of gravitinos can cause problems for the light 
element abundances of nucleosynthesis. 

On the other hand, such an increase in the binding energy would allow helium 
to form at higher temperatures, leading to an increase in the primordial helium 
abundance. This effect, neglected in PI and PI. would be superimposed on 
other effects due to changing couplings, i.e. a change in freezeout temperature 
and neutron lifetime. One can make a rough estimate of the sensitivity of the 
helium abundance to the deuteron binding energy in the following way. The 
ratio of number densities of the neutron to the proton n/p a.t the time of helium 
formation determines the primordial abundance of helium, since 99.99% of the 
neutrons go on to form helium. This is related to the initial ratio of neutron to 
proton number density at the time of weak interaction freezeout no/po, via the 
decay of neutrons into protons, by 

n Q-i-HelT 

P Pq/hq + I - 


where tne is the time at which helium production takes place and r is the neutron 
lifetime (currently measured to be 887 s). In this formula, weak interactions are 
taken to freeze out at t = 0. We have also assumed that all helium production 
occurs instantaneously at T^e- Although still a subject of some controversy, the 
orthodox view of light element abundances determines the initial mass fraction 
to an accuracy of ±5% is the ratio n/p jll]. Translating this into a constraint 
upon the variation of The we obtain 

AT' 

- 0.20 <- — < +0.13 (46) 

The 

It is reasonable to suggest that The is proportional to the binding energy of the 
deuteron, since only below a certain temperature set by that binding energy can 
helium formation proceed unimpeded by the photo-dissociation of the deuteron. 
This is because the reactions which form helium all rely upon deuterium as an 
intermediate building block :- 


H2(H2,n)He3(H2,p)He^ 

H2(H2,p)H3(H2,n)He^ 

H2(H^7)He^ . (47) 
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Using Eq. (I2H|1 we can find out how the parameter (3 changes for different values of 
the parameter c. In order to translate this into the binding energy Eb = (3Vsq it 
is necessary to take into account the scaling of the depth of the square well Vgq oc 

3 j 

oc m|A “2 so the result is not in terms of the dimensionless parameter c. 
Therefore we take the case where the variation in the quark Yukawa couplings is 
not significant and simply consider the variation in A. The result is 

- 0.04 < ^ < +0.04 (48) 

which relates into constraints on the gauge couplings at 1 TeV in non-SUSY 
models or 10^® GeV for GUT models of 


- 0.005 < 


«3 


< 0.005 

ITeV, no SUSY 


-0.0023 < 


as 


< 0.0023 

1016 GeV, SUSY 


(49) 


which is obviously very restrictive. In order to investigate this further it would 
be necessary to obtain a more detailed model for the deuteron binding energy 
and the various binding energies, cross sections and decay rates in nucleosyn¬ 
thesis. This effect then has to be added to the effect of varying neutron-proton 
mass difference and the effect of gauge coupling variation on weak coupling Hg. 
However the point remains that in general a variation in the gauge coupling at 
nucleosynthesis would have a very large effect upon nuclear binding energies, and 
consequently on the primordial abundance of light elements. This would create 
problems for models such as gg where a change in a of order 1% at nucleosyn¬ 
thesis is considered. 

If one decreases A, the binding energy of the deuteron increases, therefore 
helium is produced earlier when there are more neutrons and the helium abun¬ 
dance goes up. At the same time, a decrease in A decreases the neutron-proton 
mass difference, which increases the number of neutrons and therefore the helium 
abundance, so one would expect that two effects should work together, implying 
a yet more restrictive bound. 


6.2 Implications of a stable di-neutron or di-proton 

A negative change in the gauge coupling at high redshift would decrease the 
Goulomb repulsion and increase the nuclear forces in the di-proton system. The 
di-proton and di-neutron systems only become stable if A decreases to around 
one tenth of its present day value. This corresponds to a 10% decrease in the 
unihed gauge coupling for a high scale model and a 25% decrease for a TeV GUT. 
The stability of the di-proton would have a catastrophic effect upon the lifetime 
of stars as it would provide a rapid channel for hydrogen fusion jlU]. However 
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the constraints on the gange conpling variation in the matter dominated epoch 
are several orders of magnitnde too small for this to occur ng. The presence of 
the di-proton would also be disastrous for nucleosynthesis and might eliminate 
all the hydrogen in the universe. However, the large negative change in the gauge 
couplings required indicates that this is not a particularly strong constraint. 

It is not immediately obvious how dangerous the stability of the di-neutron 
would be to nucleosynthesis, since the neutrons would probably still end up in 
Helium atoms. Still, we can safely say that the variation in the gauge coupling 
required for the di-proton and di-neutron to become stable is much larger than 
the typical orders of magnitudes being considered at the moment. 

7 Conclusions 

We have shown in this paper that a 3% increase in the QCD coupling constant as 
at the GUT scale would result in the deuteron becoming unbound. The deuteron 
binding depends only on nuclear forces, so this conclusion cannot be escaped by 
considering the variation of more than one gauge coupling at once. Only negative 
variations at the level of 10% could bind the di-proton and the di-neutron. 

We have developed formalisms which enable one to calculate the variation 
in low energy parameters as a function of the variation of gauge and Yukawa 
couplings in the underlying theory: in many models one expects the variation 
in the electroweak and SUSY-breaking sectors to be the dominant effect at low 
energies, but the model-dependent nature of such effects means that no firm 
conclusions can yet be reached. 

We have also considered the effect of variations in the binding energy of the 
deuteron on the time at which helium formation occurs, and consequently on 
the helium abundance. This effect is complementary to the other effects on 
nucleosynthesis due to variation of gauge couplings, but on its own it constrains 
variation in at nucleosynthesis to within 0.25 %. 
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